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ABSTRACT 



The relative proximity of the recent gamma ray burst (GRB) 030329 resulted 

in a large gamma-ray fluence and in the brightest-ever afterglow (AG), hours 

after the burst, in the radio, optical and X-ray bands, permitting precise AG 

! measurements, sensitive tests of models and an excellent occasion to investigate 

the association of GRBs with supernovae (SNe). The Cannonball (CB) model 

^ { '. provides a good, simple and universal description of all AGs of GRBs of known 

redshift, so that it is straightforward to use it to predict what the expected SN 

signatures are. In the case of GRB 030329, 10 days after burst the AG should 

begin to reveal the lightcurve, spectrum and polarization of an underlying SN 

— akin to SN1998bw — which will peak in the NIR/optical band around day 15. 

These effects will be easily observable if indeed SN1998bw is a good "standard 

candle" for GRB-associated SNe and if the so far unknown extinction in the host 
O , , 

C<") , galaxy is not too large. 

O 
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It is important to know what it is that makes GRBs. The evidence for an association 



of long duration GRBs with supernovae (SNe) has been treated by observers with due care, 
since in only a few cases is the evidence very clear and model-independent. Moreover, the 
only overwhelmingly clear association — that of GRB 980425 with SN1998bw — does not 
fit at all in the generally accepted fireball models: the GRB is far too weak and the SN is 
annoyingly peculiar. These may be the reasons why authoritative reviews of the subject, 
even very recent ones, do not consider the SN/GRB association to be an established fact (e.g. 
Waxman 2003). To judge from recent Gamma Ray Burst Communication Network (GCN) 



Mado@phep3.technion.ac.il, arnon@physics.technion.ac.il, dar@cern.ch. 
Physics Department and Space Research Institute, Technion, Haifa 32000, Israel 

2 alvaro.derujula@cern.ch. Theory Division, CERN, CH-1211 Geneva 23, Switzerland 



- 2 - 



postings concerning GRB 030329, the observers' attitude is more positive: many GCNs 
discuss minor changes in the smooth afterglow (AG) decline as possible SN signatures, even 
much before — and a fluence level much larger — than one may expect them to occur. 

In the CB model long-duration GRBs are assumed to be associated with common core- 
collapse SNe (Dar and De Rujula 2000). In this model there is nothing special about GRB 
980425, which is describable in exactly the same way as all other GRBs of known redshift; 
nor is SN1998bw particularly peculiar: it did not emit the observed puzzling radio and X-ray 
signals, the CB responsible for the GRB did it, as discussed in tiresome detail in Dado el 
al. (2003a). It follows that, in the CB model, it makes sense to use SN1998bw as a putative 
standard candle for SNe associated with GRBs. The result of that exercise is that, in the 
CB model, there is evidence for a GRB/SN association in all cases in which the SN could in 
practice be observed, about 1/2 of the total (Dado et al. 2002a, b,c; 2003a,b,c). 

The Cannonball (CB) model provides a good, simple and universal description of all 
AGs of GRBs of known redshift, z, so that it is straightforward to use it to predict what the 
expected SN contribution to the AG should be, since one is confident on the extrapolated 
CB contribution, a "background" to the SN signal. By making a fit to the observed AGs, 
before the SN contribution became observable, we have, in three cases, succesfully predicted 
what this contribution would be (Dado et al. 2002b, c, 2003b). We endeavour to do the same 
for GRB 030329, by submitting this paper to the astro/ph archives ~2 days before the SN 
contribution ought to become discernible. The credibility of our predictions is weakened by 
the fact that the data, so far preliminary and unpublished, are insufficient to estimate the 
extinction in the host galaxy. Also if — unlike in all other cases for which there is evidence of 
an association — the SN that generated GRB 030329 turns out to be more than two times less 
intrinsically luminous than SN1998bw, it will be difficult to see its contribution to the AG 
light curves, though its presence may still significantly distort the late-time AG spectrum. 

We also show that -in the CB model- most of the "unusual" features of the AG of 
GRB 030329 occured before in other GRBs and could, in this sense, be anticipated. The 
early optical AG of GRB 030329 is, as for GRBs 990123, 021004 and 021121, a direct tracer 
of the expected circumburst density profile: in all these cases the AGs were "caught" very 
early and they provide extremely clean evidence for the expected 1/r 2 profile generated by 
a pre-SN stellar wind (Dado et al. 2003c). Both GRB 030329 and 021004 have 7-ray light 
curves dominated by two pulses or, in the CB model, two cannonballs. Both AGs display 
two wide shoulders that, in the same model, correspond to the contributions of the two 
CBs. The superimposed variations in the predicted smooth AG lightcurve, ups and downs 
of ~ 1/2 magnitude, are, as for GRBs 021004, 000301c and 970508, to be expected: they 
directly trace, like a needle on a vynil record, moderate deviations from a const ant- density 
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interstellar medium. 



1. GRB 030329 

The bright optical AG of this GRB (Vanderspek et al. 2003), which was first observed 
~ 1.25h after burst by Peterson and Price (2003) and by Torii (2003) is — at magnitude 
i?~12.5 — by far the brightest optical afterglow observed hours after burst. The X-ray AG, 
also very bright, was seen at time t ~ 5h after burst with RXTE, by Marshall and Swank 
(2003). The bright radio AG was observed at t~14h by Berger et al. (2003). The redshift 
of the host galaxy was first determined by Greiner et al. (2003) to be z = 0.1685. The large 
luminosity of this AG has triggered a large interest, reflected in a copious and enthousiastic 
release of GCN communications. 



The Cannonball Model of GRBs 

In the CB model (Dar and De Rujula 2000, 2001; Dado et al. 2002a; 2003a) reviewed 
in De Rujula 2002 and Dar 1993), long duration GRBs and their AGs are produced in core- 
collapse supernovae akin to SN1998bw by the ejection of bipolar jets of ordinary-matter, 
hydrogenic plasma clouds or "cannonballs" (CBs) with high Lorentz factors (70 ~ 10 3 ). A 
CB is emitted, as observed in /i-quasars, when part of an accretion disk falls abruptly onto 
the newly-born compact central object. Crossing the circumburst shells with a large 70, the 
surface of a CB is collisionally heated to keV temperatures and the thermal radiation it 
emits as it reaches the transparent outskirts of the shells — boosted and collimated by the 
CB's motion — is a single 7-ray pulse in a GRB. The cadence of pulses reflects the chaotic 
accretion and is not predictable, but the individual-pulse temporal and spectral properties 
are (Dar and De Rujula 2001; Dar 2003). In practice GRBs are observable only if the 
angle 9 subtended by the CBs' velocity vector and the line of sight to the observer is small: 
= O(l/7o). 

After becoming visible, a CB first cools by bremsstrahlung and expansion, emitting a 
hard spectrum that is seen in the X-ray band with a fluence decreasing with the expected 1/t 5 
behaviour. When a CB's temperature approaches ~ 1 eV (within a few observer minutes), 
its emissivity is dominated by synchrotron emission from the electrons that penetrate in it as 
it propagates in the interstellar medium (ISM). Integrated over frequency, this synchrotron 
emissivity is proportional to the energy-deposition rate of the ISM electrons in the CB. 
These electrons are Fermi accelerated in the CB's tangled magnetic maze to a broken power- 



-4- 



law energy distribution with a "bend" energy equal to their incident energy in the CBs' rest 
frame. Their synchrotron radiation — the afterglow — is also collimated and Doppler-boosted 
by the relativistic motion of the CBs. The radiation is also redshifted by the cosmological 
expansion and attenuated on its way to an earthly observer, during its passage through the 
CB itself, the host galaxy, the intergalactic space and our own galaxy. 



2. The NIR/Optical AG in the CB model 



In the CB model, the observed AGs have three origins: the ejected CBs, the concomitant 
SN explosion, and the host galaxy. These components are usually unresolved in the measured 
"GRB afterglows", so that the corresponding light curves and spectra are the cumulative 
energy flux density: 

Fag = F C Bs + F S n + F HG . (1) 

The contribution of the host galaxy, F HG , is usually extracted from "very" late time obser- 
vations when the CB and SN contributions become negligible, or is best fitted if such data 
are not available. 

Let the energy flux density of SN1998bw at redshift = 0.0085 (Galama et al. 1998) 
be FtnJy, t\. For a similar SN placed at a redshift z: 



FsnM = 



1 + z D 2 L (z bw ) 



l + z bw D\{z) 
1 + z 



X 



bw 



V ■ 



%bv 
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(2) 



where A$n{v, z) is the attenuation along the line of sight and D^{z) is the luminosity distance 
(we use a cosmology with VLm = 0.3, Qa = 0.7 and H = 65 km/s/Mpc). 



In the rest frame of a CB the AG that an observer sees at optical frequencies is given 



by: 



F CB M 



/[ 7 (t)] 2 n p [u/u b ] 



-1/2 



(3) 



Vb y/l + [ls/lS b ](P-V ' 

where / is a normalization constant (see Dado et al. 2003a for its theoretical estimate), 7(t) 
is the Lorentz factor of the CB, p pa 2.2 is the spectral index of the radiating electrons in 
the CB, and v b is the "injection bend" frequency for an interstellar density n p : 



v b ~ 1.87 x 10 3 [ 7 (t)f 



n r 



10 3 cm 3 



1/2 



Hz. 



(4) 
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The theoretical motivation, as well as the excellent observational support for this "bend" , are 
discussed in Dado et al. 2003a. An observer in the GRB progenitor's rest system, viewing a 
CB at an angle 9, sees its radiation Doppler-boosted by a factor 5: 

5(t) = 1 2 7ft) (5) 

1 ; 7 (t)(l-/3(t)cos0) l + 02 7 (t)2 ' 1 } 

where the approximation is valid in the domain of interest for GRBs: large 7 and small 
9. The cannonballs' AG spectral energy density Fq^ seen by a cosmological observer at a 
rcdstiift z is* 

^,,,A(i,,t)(l + z)/,{t)' \(\+z)v S(t)f 



S(t) 'l + z 



(6) 



where A(v,t) is the correction for the total extinction of the CB's radiation. 



For an interstellar medium of constant baryon density n p , the Lorentz factor 7 (i) is 
given by: 



7 = 7(7o,0, zoo; *) = B- 1 [9 2 + C9 4 + l/C] 
2/ (b 2 + 2 9 6 + B VB 2 + 4 9 & y "" 



C = 



B = l/^ + 39 2 / l0 + Qct/[(l + z)x oo ] (7) 

where 70 = 7(0), and = N C b/ Rmax n p) characterizes the CB's slow-down in terms of 
Ncb- its baryon number, and R ma x- its asymptotic radius (it takes a distance £00/70 for the 
CB to half its original Lorentz factor). 

For AGs that are observed very early (those of the GRBs 021211, 990123, 021004 and 
030329) the approximation of a constant ISM density, adopted above, is not good, since 
the CBs are first observed while travelling in the progenitor star's wind. The total wind 
grammage is insufficient to make 7 (i) deviate from 7o during the CB's traversal of the wind. 
Thus the only modification is that the density occuring in Eqs. (3) and (4) now becomes, 
for a typical "windy" profile: 

n(r) =n p (l+r 2 /r 2 ), (8) 

where n p is again a constant. These fits have two extra parameters (n p , which was previously 
embedded with /, or played a very marginal role via Eq. (4), and tq). Note that at very 
early time, r fa ct~f 5 /(l + z) oc t and the above circumburst wind profile yields F u oc 
[1 + (t/t) 2 ] 0,75 (Dado et al. 2003c), where t is the observer time when r = r ). This very 
simple expression fits like the proverbial glove the data on the early time behaviour of the 
AGs of the GRBs 021211, 990123, 021004 and 030329, that in the fireball models would be 
attributed to a "reverse" shock. 
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The contribution to A(u, z) in Eq. (6) from selective extinction in the host galaxy and in 
the intergalactic medium can be estimated from the difference between the observed spectral 
index at very early time when the CBs are still near the SN and that expected in the absence 
of extinction. Indeed, the CB model predicts — and the data confirm with precision — the 
gradual evolution of the effective optical spectral index towards the constant value — 1.1 
observed in all "late" AGs (Dado et al. 2002a, 2003a). The "late" index is independent of 
the attenuation in the host galaxy, since at t > 1 (observer's) days after the explosion, the 
CBs are typically already moving in the low column density, optically transparent halo of 
the host galaxy. 

The comparison in Dado et al. (2002a,b,c; 2003a,b,c) of the predictions of Eq. (6) with 
the observations of optical, X-ray and radio light-curves and spectra for all GRBs of known 
redshift are very simple, satisfactory and parameter-thrifty. 

3. GRB 030329 in the CB model 

Six properties of AGs in the CB model are particularly relevant to GRB 030329: 

• The optical AG is "caught" at very early times, when the CBs are still crossing the 
parent star's wind-generated density profile. 

• In high-precission data, density inhomogeneities along the CBs' trajectory — expected 
within star formation regions and upon exit from the superbubbles where most SNe 
take place — result in observable achromatic "bumps" in the AG. These features have 
been seen in GRBs 970508, 000301c and 021004 (Dado et al. 2002a, 2003b). 

• Except for GRB 021004, the past AG data were course enough or started late enough 
for the contributions of different CBs (which can often be resolved as individual pulses 
in the GRB phase) to coalesce into an AG describable by a single CB or a collection 
of similar ones. But individual CBs may have different properties, or be emitted at 
somewhat different angles, as observed in the /i-quasar SS 433 (Margon 1984). 

• The predicted AG spectra and their evolution, in particular the steepening at the time- 
varying frequency Ub(t) of Eq.(4) towards i/ -11 at late time, is very well supported by 
the data (Dado et al. 2002a, 2003a). 

• The excess polarization of AGs above that induced by the ISM in the Galaxy may be 
largely due to the host galaxy's ISM. In that case, it should be correlated with the 
extinction in the host and decline with time as the CBs move into its halo. 



-7- 



• The 7-ray light curve of GRB 030329, like that of GRB 021004, shows two prominent 
pulses (Vanderspek et al. 2003; http:// space.mit.edu/ HETE/ Bursts /GRB030329/). 

In the CB model the two 7-ray pulses of the last item correspond to two dominant 
CBs. The good quality and early start of the optical data for GRB 030329 implies that the 
individual contributions of the two CBs are discernable in the AG: we must fit the broad- 
band AG light curves with the additive contributions of two CBs, emitted almost in the 
same direction, but with otherwise free parameters (normalization, 70 and £00) ■ We fix the 
spectral index p in Eq. (3) to the theoretically expected p = 2.2, assume the density profile 
given by Eq. (8), and fit simultaneously all the reported NIR, optical and radio data. 

The individual CBs' contribution are given by Eq. (6) which implicitely uses all equations 
from to (3) to (7). The free parameters are the radius tq at which the wind and ISM densities 
are equal and, for each CB, the viewing angle, 6, the normalization constant / x n p , the 
initial Lorentz factor, 7(0) and the decelaration parameter, x^. 

The i?-band magnitude of a SN akin to 1998bw, displaced to z — 0.1685, is i?~20.2 
at peak brightness. Thus, in our fits, we neglect the smaller contribution of the host galaxy, 
i?>23.1 at 2a confidence level (Blake and Bloom, 2003; see also Wood-Vasey et al. 2003). 

3.1. NIR-Optical AG 

In Fig. 1 we show the CB-model's fit to data for the NIR-optical light-curves reported 
in the GCN circulars 1990, 2001, 2024, 2046, 2051, 2054, 2079 (Burenin et al); 1991, 2005, 
2028 (Rumyantsev et al.); 1992 (Andersen et al.); 1995 (Rykoff et al.); 1999 (Gal- Yam et al); 
2002, 2035, 2091 (Lipunov et al.); 2012 (Martini et al.); 2016 (Masi et al.); 2019 (Smith et 
al.); 2021 (Halpern et al.); 2022, 2075 (Zharikov et al); 2029 (Klose et al.); 2030 (Bartolini 
et al.); 2034, 2045, 2049, 2060 (Lipkin et al.); 2036 (Garnavich et al.); 2037 (Moran et al.); 
2040, 2096 (Lamb et al.); 2041 (Stanek et al.); 2042 (Schaefer. et al.); 2048 (Zeh et al.) 
2056, 2065, 2070 (Fitzgerald & Orosz); 2063 (Weidong Li et al.); 2066 (Tober et al.); 2050, 
2067, 2083, 2097 (Pavlenko et al.); 2068 (Bloom et al); 2074 (Cantiello et al.); 2077, 2084 
(Ibrahimov et al.); 2080 (Sato et al.); 2094 (Khamitov et al.); 2096 (Lee et al.); and 2098 
(Ibrahimov et al.). 

All the observational data was recalibrated to the field photometry of Henden et al. (2003), 
either by the observers or by ourselves. We have corrected for selective extinction in our 
galaxy (Schlegel et al. 1998): E(B-V) = 0.025 in the direction of GRB 030329. We did not 
correct for extinction in the host galaxy; the preliminary data in the GCN does not provide 
clear evidence for it. 
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The broad-band fitted parameters of the two CBs are 0[1] = 2.00 mrad, a nearly identical 
0[2] = f .95 mrad for the second CB's direction, 70 [1] = 1477, 70 [2] = 976 (implying 6 [1] = 
306, 5 [2] = 423), x^l] = 476 kpc, and Xoo[2] = 36 kpc. The density and "wind" parameters 
of Eq. (8) are n p = 0.86 x 10~ 2 cm -3 and r = 26 pc, which yield pr 2 = 9.3 x 10 13 g cm -1 for 
the progenitor-wind's grammage, similar to the canonical value for the winds of the parent 
massive stars of core-collapse SNe (Dado et al. 2003c and references therein). 

To demonstrate the real quality of the fit, we have blown up the R-band results in Fig. 2. 
In the region between t ~ 1 and ~ 5 days, the data "wiggles" by as much as 30% around 
the smoother theoretical curve. It would be easy to correct for this by assuming similar 
deviations of the ISM density input in Eqs. (3) and (8), relative to a constant n p , clearly a 
moot exercise. 

A SN1998bw-like contribution, as can be seen in Figs. 1 and 2, will be observable. Since 
the contribution of the second CB is still considerable near the peak brightness of the SN (at 
t ~ 15 days) the SN contribution appears as a shoulder in the light curves. If this contribution 
is much fainter — due to extinction in the host galaxy or a possible deviation from a so-far- 
succesful SN1998bw "standard candle" ansatz — its presence could still be established by the 
change of colours from the broad band F v ~ u~ 1A behaviour of the CBs' emission towards 
F„ ~ z/- 3 - 5±L0 , the SN spectrum. 



3.2. Optical Polarization 

In the CB model, after subtraction of the Galaxy's contribution, the polarization of 
the optical AG is mainly due to the effect of the host's ISM. It is therefore correlated with 
the extinction in the host, has a fixed angular position, and diminishes with time as the 
CBs escape into the galactic halo. Consequently near peak brightness of the associated SN 
(around 15 days after burst) the polarization of the AG should be that of the underlying 
SN. 



3.3. The GRB proper 

In the CB model, due to Doppler boosting and relativistic collimation (Dar and De 
Riijula 2000), the 7-ray fluence of a GRB viewed at a small 9 is amplified by a huge factor 
51: 



GRB 4 7r £)2 ^7> W 
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relative to E 7 , the total energy in photons emitted by the CBs in their rest system. The 
total "equivalent spherical", or would-be isotropic energy, E lso , inferred from the observed 
fluence, is a factor (<5 ) 3 larger than E 1 . In Dado et al. 2002a we deduced that the _E 7 values 
of the GRBs of known z span the surprisingly narrow 3 interval 10 44±0 3 erg, the spread in 
F GRB being mainly due to the spread in their values of 5o (deduced from the fits to their 
AGs). For GRB 030329 the CB-model fit to its broad band AG yields 5 [1] 3 = 2.86 x 10 7 
and 5o[2} 3 = 7.56 x 10 7 . The CB-model expectation from the fitted AG is E iso pa (S [l} 3 + 
<5 [1] 3 ) Ej pa 1.04 x io 52±0 - 3 erg, in agreement with the observed E lso pa 1.1 x 10 52 erg, deduced 
from its measured redshift, z = 0.168, and fluence in the 30-400 keV band ~ 10 -4 erg cm -2 
(Vanderspek et al. 2003). 

At the very early stage of 7-ray emission, CBs are still dense enough for the ambient 
protons to interact hadronically with those of the CB's ordinary- matter plasma. In the 
simplest "surface model" of the generation of the GRB proper, the energy deposited in the 
CBs by these interactions is reemitted with a quasi-thermal spectrum with a predictable time- 
dependent temperature, leading to the correct predictions of the very stable characteristic 
energy of GRBs, and of the observed spectral indices of their "Band" spectrum (Dar and De 
Riijula 2000, Dar 2003). 

The light curve of a GRB's single pulse rises abruptly as the column density of the 
material in front of a CB descends to 0(1) attenuation lengths, and diminishes thereafter 
as the density of the material the CB encounters diminishes with time. For a windy 1/r 2 
density profile of this density, the shape of a pulse at fixed energy is analytical and simple: 



with t the observer's time and Ax the proper distance travelled by the CB from the onset 
of the pulse to the position at which one 7-ray absorption length is still in front of it. In 
the rough approximation of an energy-independent 7-ray attenuation in the circumburst 
material, the shape of an energy-integrated pulse is the same as in Eq. (10). 

We have constructed a rough description of the 7-ray light curve of GRB 030329 (Van- 
derspek et al. 2003; http:// space.mit.edu/ HETE/ Bursts /GRB030329/) by summing 
two pulses, both with the shape of Eq. (10), separated by 12 observer's seconds, and with 
At[l] =6s, At[2] =4s. The result is shown in Fig. 3. The description of the GRB is quite 



F CB (t) oc 




At = 



Ax (1 + z) 
cy 5 



(10) 



3 GRBs in the CB model are much better standard candles than in the standard model (Frail et al. 2001). 
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satisfactory, even though we would have expected At [1] /At [2] = 70 [2] 5q[2]/ (70 [1] 5o[l])~ 1.1, 
as opposed to 1.5. It would be easy to find excuses to exit from this minor problem (e.g. 
we have not quoted errors in our parameters, that are made somewhat meaningless by the 
slightly erratic behaviour of the AG light curves). 

4. Conclusions 

We have shown that the data on GRB 030329 and on its early NIR-optical AG are well 
described by the CB model. This is also the case for all other GRBs of known redshift, 
including GRB 980425, whose associated SN we use as a putative standard candle. 

In the CB model, we contend, the association between long-duration GRBs and SNe 
is extremely well established: in all GRBs of known redshift (all those with z < 1.12) a 
SN1998bw-like supernova could be seen, and with various degrees of significance, it was 
seen. SN1998 stood the "standard-candle" test with surprising precision. In many cases, 
though, the test depended on very significant extinction corrections in the host galaxy. If this 
extinction turns out to be significant for GRB 030329, the direct observation of a SN signal 
in the AG could be jeopardized. In such a case one would need precise colour photometry 
and spectroscopy with the most powerful telescopes, such as VLT, Magellan, Keck and HST, 
and very late time measurements of the host galaxy's magnitude. 
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Fig. 1 — The NIR-optical observations of the AG of GRB 030329 and the fit for two CBs 
with different parameters, as described in the text. The ISM density is a constant plus a 
"wind" contribution decreasing as 1/r 2 . The various bands are scaled for presentation. The 
data was selected from the GCN notices quoted in the text, recalibrated with the observations 
of Henden et al. 2003. The host-galaxy's contribution was neglected. The individual bands 
have been rescaled for clarity. 
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Fig. 2. — Blow-up of the R-band results of Fig. 1. The ISM density was assumed to be a 
constant plus an additional "wind" contribution decreasing as 1/r 2 . The wind contribution 
is only significant at t < 0.1 days, after which the CBs are more than 10 pc away from the 
progenitor. This "wind" is also seen in other AGs observed early enough (Dado et al. 2003c). 
The individual contributions of the two CBs and of a SN akin to SN1998bw (at the GRB's 
redshift) are also shown. 
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Fig. 3. — The 7-ray light curve of GRB 030329 (the red binned curve), and its simple 
CB-model description (the blue continuous line). 



